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Abstract. Thermal modeling of the LISA gravitational reference sensor (GRS) includes such ef- 
fects as outgassing from the proof mass and its housing and the radiometer effect. Experimental data 
in conditions emulating the LISA GRS are required to confidently predict the GRS performance. 
Outgassing and the radiometer effect are similar in characteristics and are difficult to decouple ex- 
perimentally. 

The design of our torsion balance allows us to investigate differential radiation pressure, the 
radiometer effect, and outgassing on closely separated conducting surfaces with high sensitivity. A 
thermally controlled split copper plate is brought near a freely hanging plate-torsion pendulum. We 
have varied the temperature on each half of the copper plate and have measured the resulting forces 
on the pendulum. 

We have determined that to first order the current GRS model for the radiometer effect, out- 
gassing, and radiation pressure are mostly consistent with our torsion balance measurements and 
therefore these thermal effects do not appear to be a large hindrance to the LISA noise budget. 
However, there remain discrepancies between the predicted dependence of these effects on the tem- 
perature of our apparatus. 
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INTRODUCTION 

The low frequency end of the LISA gravitational wave sensitivity requires that the 
proof masses are kept in free-fall below 8a = 3 x 10~ 15 (m/s 2 )/v / Hz from 0.1 to 
lOmHzllLlalEl. The spurious accelerations linked to random thermal fluctuations have 
been allocated one-third of this total budget for meeting the desired LISA sensitivity 
[4]. Thermal fluctuations on the electrode housing of the GRS are estimated to be 
Sat = 10 -10 K 2 /Hz B2D, which combined with theoretical understanding of the thermal 
effects yields baj = 3.5 x 10~ 16 (m/s )/a/Hz, which is one-third of the allocation. 
The margin in the budget that this assessment creates could be used to relax other 
requirements. Therefore it is imperative that the models which were used to generate 
the current assessment are validated in a laboratory setting. 

The three largest thermal effects in the LISA noise budget are radiation pressure f the 
radiometer effect, and asymmetric outgassing. These three terms add coherently: pLS] 
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where o is the Stefan-Boltzmann constant, c is the speed of light, M = 1.96kg is the 
mass of the LISA proof mass, A = 2.12 x 10~ 3 m 2 is the area of one face of the LISA 
proof mass, P = 10~ 5 Pa is the residual gas pressure at the proof mass, and T = 293 K is 
the temperature. The first term in Equation \T\ is due to fluctuating asymmetric radiation 
pressure, 5<3rp. The second term is due to the radiometer effect, Scire- The pressures on 
opposite faces of the proof mass equilibrate quickly so that there is nearly no pressure 
difference across it. This implies that the contribution due to the radiometer effect should 
be reduced |0]. Leading to a pressure reduction factor of dp ~ 0.1 [8]. The last term 
is due to fluctuations in asymmetric outgassing on opposite sides of the proof mass. 
The outgassing rate Q{T) most likely will follow an exponential law in temperature 
referenced to an activation temperature 0. The value of this activation temperature 
should be no more than 30000 K [3]. The value for the nominal outgassing rate is taken 
to be Qq « 4 x 10 10 Torr L cm" 2 s" 1 » 5 x 10~ 7 kg/s 3 H. The factor, G, is related to 
gas particle conductances about the proof mass and through holes in the housing walls, 
and to the outgassing area of the proof mass. The nominal value for LISA is G = 0.4m/s 
S]. 

The derivations of these three effects were made assuming a constant temperature 
gradient across the housing which is then varied. This implies that the three effects 
described above are dependent on the frequency of fluctuations only through the level of 
the thermal fluctuations themselves, v^at- 

Our apparatus provides the unique opportunity to investigate the gap spacing depen- 
dence as well as the pressure and linearity of the thermal effects. Our results are com- 
plimentary to those that have been presented by the University of Trento group [@, @]. 
Whereas they have designed a mock-up of the LISA proof mass and electrode housing 
and therefore emulate the geometry of LISA, we have an increased torque sensitivity, the 
ability to change the gap spacing, and a simpler geometry which allows ease in modeling 

iE3]. 



TORSION BALANCE APPARATUS 



Our torsion balance apparatus is described in [10]. Each of the two halves of the copper 
plate has embedded a heating element and a temperature sensor. Further temperature 
sensors are located on the plate translation stage, near the control electrodes, on the 
autocollimator, on the outside of our vacuum chamber, on the inside of the styrofoam 
enclosure, and the temperature controlled laboratory. 

We predict values for the three thermal noise sources under study for our torsion bal- 
ance tests. Since our torsion pendulum measures torques we must convert the measured 
torques into LISA-equivalent accelerations. As mentioned in ||10(l . this conversion is 
done by using a moment arm at which the spurious forces act on the pendulum and the 
mass of one LISA proof mass. The base pressure of our experiment is ps 1.2 x 10 Pa. 
We assume that the outgassing activation temperature and the outgassing rate of copper 
is similar to that assumed for the LISA GRS. The cross-sectional area overlap of one half 
of the copper plate and our pendulum is Auw = 4.7 x 10~ 3 m 2 . The gas conductance rate 
for our setup should be similar in value to those for the LISA GRS, likewise we assume 
the same value for the factor G. Note that simple gas dynamics implies that the conduc- 



TABLE 1. Expected Values of Thermal Noise Sources for LISA using the values given in 
the introduction, and LISA-equivalent accelerations for our apparatus using the values given in the 
text. The mean temperature of the copper plate and the residual gas pressure differentiates each 
column. The low and high pressure data were computed with base pressures of 1.2 x 10~ 5 Pa and 
1.6 x 10~ 4 Pa respectively. Values in the table are in units of 10~ n (m/s 2 /K) \/Sat. 

Noise Source LISA UW (25°C) UW (28°C) UW (25°C) UW (28°C) 

low pressure high pressure 
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tance rate should increase if the separation between the pendulum and plate grows, and 
therefore outgassing effects should decrease with plate-pendulum separation. 

In Table Q] we list the predicted values for the three thermal effects in Equation (OQ) 
for LISA and the predicted values for our experiment. In this table we have computed 
our LISA-equivalent noise levels when the mean temperature of the copper plate in 
our apparatus is 298 K and when it is 301 K. From these values we predict a thermally 
induced acceleration noise of 3 x 10~ n (m/s 2 /K) v^at relatively independent of the 
temperature of the copper plate. In addition, we show values for the system at an elevated 
pressure of 1.6 x 10 4 Pa, for which we expect an increase by nearly a factor of 2. 

We present our results in terms of a thermal-to-torque transfer function. The conver- 
sion from torque to LISA-equivalent acceleration is ~ 15(m/s )/(Nm), such that the 
expected thermal-to-torque transfer function should be 2 x 10~ 12 Nm/K when the sys- 
tem is at low pressure, and about 4 x 10 Nm/K for the high pressure scenario. 

The copper plate reached 35° during the two day mild bakeout of our apparatus under 
vacuum. Our measurements were made 150 days after this bakeout occurred. The copper 
plate is made of oxygen-free high-conductivity (OFHC) copper, which has an outgassing 
rate 10 times lower than that assumed in the above calculations 

THERMAL MEASUREMENTS 

The temperature difference between the two halves of the copper plate emulates a 
temperature gradient across the LISA proof mass. Our procedure is to vary this gradient 
recording the torque, thereby the LISA-equivalent acceleration, as a function of a variety 
of adjustable variables, such as the temperature gradient, the mean temperature of the 
copper plate, the plate-pendulum separation, and the residual gas pressure inside our 
apparatus. 

DC Measurements 

Applying constant power to the heater in either the left or right section of the copper 
plate shows a linear rise in temperature of that half as well as a linear response of 




FIGURE 1. Heating of one half of the copper plate (left) followed by heating of the other half some 
time later (right) in equals steps of heater output power. The top chart shows the temperature of each half 
of the copper plate as well as the mean temperature all in °C, the middle chart shows the temperature 
gradient between the two halves in K, and the lower chart shows the measured torque on the pendulum. 
The plate-pendulum separation for this data was 1 mm. 



the torque on the pendulum. Figure Q] shows heating of half of the copper plate in 
constant steps of heater power. This data was taken with a plate-pendulum separation 
of 1 mm. The heated half changes by about 0.76 K over the course of this data set. 
Thermal coupling from one half of the copper plate to the other is ~30% such that 
the other half gained 0.23 K. The temperature sensor placed on the electrode side of 
the apparatus (inside the vacuum chamber) rose by 47 mK over the course of this data 
set. Temperature sensors attached to the outside of the vacuum chamber changed in 
temperature by ~46mK. The temperature of the styrofoam enclosure did not change 
appreciably. This indicates that the thermal feedthrough from the copper plate to the 
apparatus is on the order of 10%. 

It can be seen from the slower rise in the torque as compared with the temperature 
gradient that there is a delay between the temperature gradient and the reaction of 
the feedback loop of ~ 1700 s. We suspect that this delay may be due to an effective 
activation time of outgassing: the trapped gas at the surface of the copper immediately 
escapes whereas the trapped gas within the copper slowly diffuses. Additionally, transfer 
of heat to the pendulum should raise the pendulum temperature by a fraction of the plate 
temperature change. The heating of the pendulum should increase outgassing from the 
silicon. There is a delay between the heating of the copper plate and the heating of the 
pendulum leading to a delayed response in the feedback torque. 



The temperature gradient between the two halves of the split copper plate is the 
important quantity under investigation. The change in temperature gradient makes a 
change in all three thermal sources under investigation. This causes a pronounced change 
in the torque. The thermal-to-torque transfer function for the data on the left in Figured] 
is about 2.3 pNm/K. 

After some time we heated the other half of the copper plate and brought the temper- 
ature gradient back to nearly zero. This data is on the right in Figure \T\ The thermal- 
to-torque transfer function for this data is 3.0 pNm/K. Notice that the mean temperature 
of the copper plate in the data on right is ~1.5 K higher than the data on the left. This 
increase in the transfer function appears to be a function of the mean temperature of the 
copper plate as described in more detail below. 



AC Measurements 

The important frequency band of thermal measurements for LISA is from 0.1 to 
1 mHz. Looking at coherent signals allows us to determine the thermal frequency re- 
sponse of our system. By sweeping in frequency from 0. 1 to 1 mHz we can map the 
transfer function from thermal input to torque output. Figure [2] contains data for such a 
frequency sweep at two mean plate temperatures and two gas pressures. This data was 
taken with a plate-pendulum separation of 1 mm. 

The thermal mass of the copper plate creates a low pass filter so that the temperature 
gradient amplitude is smaller at higher frequencies. Our typical amplitude is nearly 
200 mK at 0.1 mHz. From a series of measurements we found that the thermal-to- 



torque response is independent of the signal amplitude at any particular frequency 111211 . 
However, we still observe an apparent low pass response with a corner frequency near 
0.1 mHz and a slope index of —0.75 in amplitude. 

From our DC measurements we found a time delay between the temperature gradient 
and the control torque of about 1700 s. At 10 /iHz this would be a phase delay of 
— 18° which is in good agreement with the data presented in Figure [2l This phase 
delay becomes larger for frequencies above 1 mHz following the decrease in transfer 
function. However, it appears to turn around at higher frequencies. We suspect this is an 
experimental artifact due to electrical coupling in our system. 

To verify that this frequency and phase dependence was not dependent on our fre- 
quency sweeping function, we used a variety of different sweep functions, such as linear, 
logarithmic, and more complicated functions of time, as well as monochromatic sinu- 
soids. We found no dependence on the sweep function. The individual data points in 
Figure [2] are for fixed-frequency sine waves. 



DISCUSSION 

When the mean temperature of the copper plate is elevated by a few degrees we see an 
increase in the thermal-to-torque transfer function, which is not predicted by TableCD An 
increase in the outgassing rate is likely the cause. As mentioned above, it is also possible 
that the pendulum is appreciably changing in temperature such that the outgassing rate 
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FIGURE 2. Thermal-to-torque transfer function and phase delay for thermal frequency sweeps between 
30/iHz and 1 mHz. The solid (dashed) transfer function traces were taken at a base pressure of 1.2 x 
10~ 5 Pa (1.6 x 10~ 4 Pa). Raising the gas pressure increased the transfer function by a factor of 13. Raising 
the mean temperature of the copper plate by 3 K raised the transfer function by a factor of 1 .5 irregardless 
of the gas pressure in our apparatus. 



from the silicon increases as well. The outgassing rate should decay with time as a power 
law II 1 lL 1 1 3M . Observing this decay would yield support for this hypothesis. 

The radiometer effect is linearly proportional to the gas pressure in the system. 

Figure [2] contains thermal-to-torque traces at an elevated pressure of 1.6 x l<T 4 Pa. 
The thermal-to-torque transfer function at frequencies above 0. 1 mHz appears to be a 
factor of 13 higher when the system is at this higher pressure. This linear response is 
exactly expected for the radiometer effect, if the radiometer effect accounted for the 
bulk of the thermal-to-torque transfer. As demonstrated in Table [H the radiometer effect 
should be a small contribution to the total. The three contributing thermal effects add 
linearly since they are coherent in the temperature gradient. If outgassing and radiation 
pressure are independent of gas pressure, then the increase in the transfer function when 
the system is at high pressure would be due to the radiometer effect. This would imply 
that the radiometer effect is accounting for the bulk of the transfer, and therefore is a 
factor of 15 higher than anticipated, i.e., the pressure reduction factor, dp, should be an 
amplification of 1.5 not a reduction of 10. However, the radiometer effect is inversely 
proportional to the mean temperature and therefore does not account for the increase of 
1.5 when the copper plate is at an elevated temperature. 

If instead outgassing is the cause of the increase when the system is warm, then the 
outgassing contribution must also increase when the system is at an elevated pressure 
which is currently not included in the LISA models Q5Q. 



DIRECTIONS FOR FURTHER STUDY 



Our results are reasonably consistent with the current theoretical understanding of the 
thermal effects used in modeling for the LISA GRS. However, there does remain an 
anomalous increase in the thermal-to-torque transfer function when the copper plate is 
heated by a small amount regardless of base pressure. We are currently conducting mea- 
surements to help explain this. In addition, we are in the process of making measure- 
ments over a range of plate-pendulum separations to determine the distance dependence 
of the thermal effects. 
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